Efficient on-chip source of microwave photon pairs in superconducting circuit QED 
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We describe a scheme for the efficient generation of microwave photon pairs by parametric down- 
conversion in a superconducting transmission line resonator coupled to a Cooper pair box serving 
as an artificial atom. By properly tuning the first three levels with respect to the cavity modes, 
the down-conversion probability may become higher than in the most efficient schemes for optical 
photons. We show this by numerically simulating the dissipative quantum dynamics of the coupled 
cavity-box system and discuss the effects of dephasing and relaxation. The setup analyzed here 
might form the basis for a future on-chip source of entangled microwave photons. 



Introduction. - The generation of photon pairs by para- 
metric down conversion (PDC)[lJ, |3, Is |5j represents 
one of the basic ways to create nonclassical states of the 
electromagnetic field, which has found numerous appli- 
cations so far. The conditional detection of one of the 
photons enables the production of single photon Fock 
states El Furthermore, PDC is the primary method 
to generate entangled pairs of particles. Apart from the 
possibility of testing Bell inequalities 0,0jIEI3j this rep- 
resents a crucial ingredient for a multitude of applications 
in the field of quantum information science, ranging from 
quantum tclcportation through quantum dense coding to 
quantum key distribution 

With the advent of superconducting circuit quantum 
electrodynamics 0, ^J, it will now be possible to take 
over many of the concepts that have been successful in 
the field of quantum/atom optics and transfer them to 
the domain of microwave photons guided along coplanar 
waveguides on a chip, interacting with superconducting 
qubits [ill H2I IT3L Il4| ]. Recent experiments have real- 
ized the strong-coupling limit of the Jaynes-Cummings 
model known in atom optics, employing a superconduct- 
ing qubit as an artificial two-level atom, and coupling it 
resonantly to a harmonic oscillator (i.e. a cavity mode[l(l| 
or a SQUID 0]). Dispersive QND measurements of the 
qubit state, Rabi oscillations and Ramsey fringes have 
been demonstrated [161 Il7| , leading to a fairly detailed 
quantitative understanding of the system, which behaves 
almost ideally as predicted by theory 

In this paper, we will analyze a scheme that imple- 
ments parametric down-conversion of microwave pho- 
tons entering a transmission line resonator coupled to 
a Cooper pair box (CPB) providing the required non- 
linearity (Fig. 0. This represents the limit of a single 
artificial atom taking the place of the nonlinear crystal 
usually employed in optical PDC experiments 
with the cavity enhancing the PDC rate (cf. 
contrast to other solid state PDC proposals [19|, 
both the basic cavity setup and the possibility of eject- 
ing the generated photons into single-mode transmission 
lines with a high degree of reliability are already an exper- 
imentally proven realitv [l(l llfil Il7j . Recently, squeezing 




Figure 1: (color online) Schematic setup for the proposed 
parametric-down conversion (PDC) experiment in supercon- 
ducting circuit cavity electrodynamics, with a Cooper-pair 
box (CPB) interacting with the three lowest modes of a trans- 
mission line resonator, whose voltage distributions are shown. 



and degenerate parametric down-conversion have been 
analyzed theoretically for a circuit QED setup cou- 
pling a charge qubit to two cavity modes. However, un- 
like the experiments and most of the theoretical investi- 
gations mentioned above, in this paper we propose to go 
beyond the regime where the box may be regarded as a 
two- level system (qubit), making use of its first three lev- 
els. By further employing its advantage over real atoms, 
namely its tunability via the applied magnetic flux and 
the gate voltage, this enables us to bring the transitions 
between the first three box levels into (near) resonance 
with the first three cavity modes (Fig. GJ, thereby drasti- 
cally enhancing the resulting probability of (nondegener- 
ate) parametric down-conversion |3w) 1— * \uj) \2lu). This 
represents the major advantage of the present scheme. 
We treat the full quantum dissipative dynamics of the 
box-cavity system, incorporating the radiation of pho- 
tons from the cavity as well as nonradiative decay pro- 
cesses and dephasing in the CPB. We will present results 
for the down-conversion efficiency, discuss the minimiza- 
tion of unwanted loss processes, and comment on possible 
applications in the end. 

The model. - The CPB is a device 0] in which Cooper 
pairs can tunnel between two superconducting islands 
due to a Josephson-coupling Ej (tunable by an external 
magnetic flux in a split-junction geometry). The number 
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Figure 2: (color online) Left: Energy- level diagram of the 
CPB - transition energies in units of Ec, for Ej/Ec = 3, as 
a function of gate charge Nq- At particular gate values (away 
from the "qubit" regime near No = 1/2), the CPB transitions 
frequencies are related in an integer ratio, (e2 — = 1:1 

or 2 : 1, respectively. This enables to match them with the 
cavity modes, giving rise to particularly efficient degenerate 
or non-degenerate PDC (2u Hiu + worSunu-l- 2ui, re- 
spectively). Right: Simplified transition scheme for the non- 
degenerate PDC process considered in the text, with a detun- 
ing A of the intermediate state |1, 2lo): e\ — ui + A. 



of transferred Cooper pairs N determines the charging 
energy, whose scale Eq = e 2 /2Cs is set by the total ca- 
pacitance of the box island, and which can be controlled 
by application of an external gate voltage (expressed in 
terms of a gate charge Nq)- 



H C pb = 4E C (N-N G ) 2 

Z N 



(1) 



Here \N) represents a charge state of the CPB. Two 
gaps in a superconducting coplanar waveguide act as 
mirrors of a cavity for microwave photons, -ff C avity = 
Sj=i WjQjQji where we will focus attention on the three 
lowest-lying cavity modes with ojj — ju> (we set h = 1). 
The electric field inside the cavity acts on the CPB via 
a capacitive coupling, adding a fluctuating quantum- 
mechanical component to the gate charge Nq. This re- 
sults |9j in an interaction 
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The coupling constants (vacuum Rabi frequencies) gj = 

goVj<t>j( x )> with go — are given in terms 

of the mode functions 4>\{x) = sin(x7r/L), 4>2{x) — 
cos(2ttx/L) and ^(x) — sn\{Zirx/ 'L), which are defined 
on the interval x = —L/2 . . . L/2. 

The full Hamiltonian forming the basis of our analysis 
is thus given by 



H = H, 



cavity 



■H, 



CPB 



H + H e 



(3) 



where H onv includes the coupling to the various reservoirs 
forming the environment. This involves both the possi- 
bility for microwave photons to leak out of the cavity, as 
well as the various possible nonradiative decay and deco- 
herence processes acting on the CPB. The details will be 
specified below when setting up the master equation. 

Basic considerations. - At least three basic features 
distinguish such a setup from the usual optical pho- 
ton PDC experiments employing nonlinear crystals: (i) 
There is no momentum conservation, as the system is 
essentially zero-dimensional, (ii) However, energy con- 
servation is much more restrictive, as the set of possible 
frequencies is limited to the discrete cavity modes. For 
appropriate input frequency, this results in a resonant en- 
hancement of the PDC process. In contrast to a passive 
filtering scheme, the bandwidth of the generated photons 
is reduced without diminishing the signal intensity, (iii) 
The microwave polarization is fixed and thus cannot be 
used for entanglement. At the end of this paper we will 
point out other options that can be explored. 

Estimating the PDC rate. - If the CPB is operated as 
a two-level system (qubit) [22j|, the decay of a 3lu pho- 
ton into two lower-energy photons comes about through 
multi-step transitions, where at least one of the interme- 
diate virtual states will have an energy detuning of the 
order of u, which contributes a small factor (g/oj) 2 to the 
PDC rate. 

However, we can enhance the PDC rate by exploiting 
more than the first two levels of the CPB. Indeed, by 
properly tuning the Josephson coupling Ej and the gate 
charge Nq, it is possible to make the transitions between 
the first three CPB energy levels |0) , |1) , \2) resonant 
with the cavity modes. The PDC process we want to 
consider is thus |0,3w) i— > \2) i— > |1, 2o>) i— » |0,w,2w), 
with £2 ~ 3w and t\ = lo + A, such that all the tran- 
sitions are (nearly) resonant. This reduces the largest 
energy denominator to the detuning A, resulting in an 
enhancement of the PDC rate by a factor of (w/A) 2 . 
What limits the enhancement? If A is made too small, 
one may end up with less than a complete pair of down- 
converted photons: Instead of passing virtually through 
the intermediate state |1, 2u) containing a 2w-photon and 
the CPB in its first excited level, that state will acquire 
a significant real population. As a result, the tempo- 
ral correlation between photons would be destroyed, and 
nonradiative decays |1, 2lj) h- > |0, 2u>) may occur, without 
emitting the second photon of frequency u>. Clearly, there 
is a tradeoff between the achieved PDC probability and 
the fidelity of down-conversion. This will be confirmed 
by the detailed simulations below. 

We will find that it is possible to achieve a PDC proba- 
bility in the percent range that surpasses that of the most 
efficient modern optical PDC schemes [i^j (which have 
a PDC probability of about 10~ 4 ; though the absolute 
PDC rate in those experiments is about 10 9 times larger 
due to the drastically higher input power). Earlier well- 
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known optical PDC experiments |5| generated less than 
one usable coincidence detection event for every 10 13 in- 
coming photons. 

Simulation of the quantum- dissipative dynamics. - In 
the ideal case, one could integrate out the intermedi- 
ate state, yielding an effective PDC term of the form 
|0) (2| d\a\ + h.c However, here we take into account 
all loss processes, by solving for the full dynamics of the 
CPB/cavity-system under an external microwave drive 
of frequency uj- m ~ 3lu, using a Markoff master equation 
of Lindblad form: 

= (^0 + £drivc + ^cavity + ^CPB + ^CPb)/ 5 ( 4 ) 

Here Cop = —i[Ho, p], and the external microwave input, 
at a frequency u)- m rs 3w and with an amplitude a, is 
described by £drivc/5 = -i [if drive (*),/>], with H drivc (t) = 
aa\e- iu)i ^ + h.c. 

The dissipative terms in the Liouvillian are of Lindblad 
form 

C[A]p = A0 (5) 

They describe: the decay of each cavity mode at a rate 
k 3 (j = 1,2,3), 

3 

pure dephasing processes in the CPB that do not lead to 
transitions between levels (at rates 7 VJ ), 

3 

and nonradiative relaxation processes leading from a level 
I to a lower energy level j of the qubit: 

^=E^«)ao (8) 

We keep only resonant terms ("rotating wave approx- 
imation", RWA) in the CPB-cavity interaction and go 
over to a frame rotating at w; n , which eliminates the 
time-dependence in ifdrive, but replaces Hq by Ho — W, 
with 3W/uj in = J2j ja)aj + |1) (1| + 3 ]2) (2|. This is 
accomplished by applying the time-dependent unitary 
transformation exp(iWt)A exp(— iWt) to the Hamilto- 
nian Ho + i?drivc and the density matrix. 

We have obtained numerical solutions of the master 
equation for a wide range of parameters, and the results 
are shown in Figs. 01 and All these simulations have 
been performed in a Hilbert space that has been trun- 
cated under the assumption of a small external drive 




Figure 3: (color online) Parametric down-conversion prob- 
ability Ppdc = P(l w ) 2w))K 2 /|a| 2 , as a function of the mi- 
crowave input frequency Wi„ and the detuning of the interme- 
diate state. The dashed line indicates the analytical resonance 
condition (see main text), the full line denotes the location of 
minimal Qi (high fidelity of PDC), and the dotted line is the 
cross-section shown in Fig. ^1 



(The maximum excitation energy of the qubit+cavity 
system is restricted to 3a; , and the down-conversion rate 
is linear in the input power). For the plots shown be- 
low, we have used experimentally reasonable parameter 
values: go/co = 10~ 2 , k 3 /uj = 10~ 4 , 7 3 - < _j/cj = 1CT 5 (j < 
I), Jip.j/u — 2 • 10~ 4 (j > 0). For reference, note that 
u ~ 10 GHz in typical experiments^]]. We have placed 
the CPB at x = 0.3 x L. 

Discussion. - In order to interpret the results, we note 
that the production of photon pairs at a rate Tprjc is 
balanced by the decay of photons out of the cavity, at 
a rate k. Thus, in an ideal lossless cavity PDC scheme, 
the probabilities to find the cavity in the states |w,2a;), 
\u>) and \2u) all become equal to Fpdc/(2k). Therefore, 
we define Fpdc/(2k) = P^^ui)- The rate R of incoming 
photons is given by R — 2\a\ 2 / 'n. Thus, the PDC prob- 
ability (chance of a given photon undergoing PDC) can 
be set to P PD c = r PDC /i? = P\ u ^u)^ 2 /M 2 (we will use 
this as a definition even where the scheme deviates from 
ideal conditions, discussing the PDC fidelity separately). 

In Fig. |3 we have plotted the PDC probability as 
a function of the input frequency Win and the detuning 
A = 6\ — lu. In general, Ppdc becomes maximal when 
Lu- m matches the doublet frequencies and 3cj (modi- 
fied by the vacuum Rabi splitting): See the two vertical 
"ridges" , independent of A. There is a third resonance 
at Win = 3w— <7i/A, with a dispersive shift depending on 
A, for which the input frequency matches the energy of 
the outgoing state (dashed curve). Here we have defined 
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Figure 4: (color online) Top: The parametric down con- 
version probability as a function of input frequency along 
the cross-section indicated in Fig. |3J and the probabili- 
ties P„ = fP(\uj)) and P 2 ^ = fP(\2u}) for one or two 
photons in the cavity, rescaled in the same manner, with 
/ = K 2 /\a\ 2 (ideally P PD c = P u = Pau). The probability 
of an excited qubit, P2 = fP(\2ui)), is also shown. Bottom: 
The "non-ideality parameters" Qi — 7-P|i,2^)/( K -P|w,2^)) an d 



Q2=P\ U) /P W 



- 1, reaching a minimum at the middle peak 



(thin line: for half the dephasing rate of 7^/aj = 2-10 



9l = ,9i 



N 



as the actual vacuum Rabi frequency. 
Although Ppdc becomes maximal when any of these res- 
onances cross, this does not yield ideal photon pairs, as 
will become clear shortly. Turning to Fig. El (top), we see 
a cross-section of Ppdc (as a function of tJi n ). In addi- 
tion, we have included the probabilities of having one or 
two photons inside the cavity, which (in this normaliza- 
tion) should be identical to Ppdc in an ideal scheme. The 
apparent surplus in w-photons is due to a nonradiative 
process, as we will discuss now. 

The relevant decay routes are: (i) |2) i— > |1) <-> |o;) i— > 
|0), (ii) \u,2u) i-> \oj) <-> |1) i-> |0), (hi) \l,2w) h-> |1) h-> 
|0), and (iv) |1, 2u>) i— > \2u) i— > |0). Process (i) leads to a 
single uj photon being emitted, while (ii)-(iv) produce a 
single 2lu photon with no corresponding partner photon. 
All of these processes get suppressed with an increasing 
energy mismatch |A| = \e\ — u)\ between the states |1) 
and \uj). A larger broadening of the levels (produced 
by dephasing or decay) may partially overcome this en- 
ergy mismatch, leading to a higher rate of unwanted loss 
processes. In order to quantify these processes, we have 
plotted, in Fig. 0] (bottom), the "non- ideality measures" 

Ql = 7-P|l,2cj}/(«:-P|a;,2 w )) and Q 2 = P\ u )/P\u,2u) ~ 1 • 

Here Qi measures the ratio of nonradiative relaxation 
from the intermediate state |l,2w) to the rate of pair 
emission, while Q2 indicates the degree to which the pop- 
ulations of \uj) and \u>, 2u) are identical (which is the case 
ideally, for kj = n). Both Qi and |Qa| should be as small 
as possible. The doublet peaks yield a large PDC rate, 
but also a large population of the CPB excited state |2), 
leading to the decay process (i) and a resulting surplus of 
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Figure 5: (color online) (a) Parametric down-conversion prob- 
ability Ppdc, and (b) the non-ideality parameter Q2, as a 
function of the detunings between the Cooper-pair box en- 
ergy levels ei,2 and the cavity modes (controlled by Ej and 
Nq). At each point, the microwave input frequency LOi n has 
been chosen to minimize the parameter Q\ (cf. Fig. 2J. The 
dotted line indicates the cross-section displayed in Fig. [|] 



w-photons (Fig. top). Thus, we observe that the vac- 
uum Rabi splitting of the doublet is essential: it allows 
for the appearance of the third (middle) peak in Ppdc 
that has a far lower qubit population and corresponding 
rate of unwanted loss processes (minima in Q\,2)- Any 
reduction in the broadening of the peaks (set by k, 7, "f v ) 
helps to further increase the quality of PDC. 

The PDC quality can also be measured by 
evaluating the 2-photon correlator Kji(t) = 

(&jaj(t)aj(t)ai\ I (J^ajai) ^a]aj-^J , which deter- 
mines the probability to detect a mode j photon at 
time t inside the cavity, provided a mode I photon has 
been detected at t = 0. Using the quantum regression 
theorem applied to our master equations, we have 
checked that at small values of Qi,2, the ideal case 
is indeed observed, where the correlator decays at a 
rate k, both for (J, I) — (1,2) and {j, I) = (2,1), while 
nonradiative processes change these decay rates and 
make them unequal. 

The gate charge Nq and the Josephson coupling Ej 
change the relevant energies ei and €2 of the first two 
excited states of the Cooper pair box, and, to a lesser 
degree, the matrix elements for the coupling to the cavity 
field. At a given value of Ej, one can tune to = 
Nq[Ej], such that the bare resonance condition €2 — 
3lj is fulfilled. In the following, we consider the effects 
of a small additional "offset gate charge" SNq, which 
mainly changes £2, while Ej itself is used to tune e\. 
The plots discussed so far have been obtained at fixed 
5Nq, while changing Ej. At any given value of Ej and 
5Nq, it is possible to select an input frequency u- m which 
minimizes Q\. The resulting PDC rate has been plotted 
as a function of the detunings of the two CPB energy 
levels, ei and £2, in Fig. 03 (a), while the corresponding 
parameter Q 2 is shown in Fig. 0(b) (Qi is below 0.01 in 
the relevant region where Q2 is small). 

We conclude that, while maintaining a good reliability 
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of the PDC process, the down conversion probability can 
become on the order of a few percent in the present setup, 
which thus indeed represents a highly efficient source of 
photon pairs. Similar results have been found for other 
reasonable parameter sets (e.g. ga/uj = ICE 2 , Kj/co = 
10- 5 , 7i ^/c; = 1CT 4 (j < I), y Vtj /u> = 10- 3 (j > 0)). 

Generation of entanglement. - The down-converted 
lu— and 2tj-photons can independently leak out of either 
side of the cavity. By post-selecting (cf. 0, onr y 
events where a photon is detected both in the left and the 
right arm each, one ends up with a frequency-entangled 
state that is directly equivalent to the entangled triplet 
state: \2u>) L Cg> \w) R + \u) L ® \2w) R .We note, however, 
that a full Bell test requires measurements in a superpo- 
sition basis, which is hard to realize for states of different 
energies. 

Another, simpler, possibility is to test for energy-time 
entanglement as first proposed by Franson 0, 0, . 
This requires feeding the generated photons into Mach- 
Zehnder interferometers, each of them containing a short 
and a long arm as well as a variable phase-shifter in one 
of the arms. By measuring the photon-detection corre- 
lation between the altogether four output ports of the 
two interferometers, it is possible to violate the usual 
kinds of Bell inequalities. The great advantage of such 
a scheme (particularly in the context of superconducting 
circuit QED) is that it does not require the polarization 
as a degree of freedom. 

A less demanding, first experimental test of the PDC 
source described here might measure the intensity cross- 
correlations of the microwave output beams (at u) and 2u> ) 
or implement homodyning techniques 6, 26] to charac- 
terize the quantum state. Finally, it is worth noting that 
for u>3 ~ €2 a sufficiently strong vacuum Rabi splitting 
between |3o;) and \2) in principle enables a scheme where 
a Rabi n pulse is used to put exactly one excitation into 
the system, which then decays in the way described here, 
thus realizing a source of microwave photon pairs on de- 
mand. 

Conclusions. - In this paper we have described and 
analyzed a setup for parametric down conversion in su- 
perconducting circuit cavity QED, suitable for the gen- 
eration of pairs of entangled microwave photons. In con- 
trast to earlier discussions, we have considered employ- 
ing a transition via the first three levels of the artificial 
atom (Cooper pair box), which can be tuned to achieve 
a drastically enhanced PDC rate. We have analyzed the 
trade-off between optimizing the PDC rate and minimiz- 
ing loss processes, by carrying out extensive numerical 
simulations of the quantum dissipative dynamics. The 
setup described here can be realized by moderate mod- 
ifications of existing experiments, and it can hopefully 
form the basis for more detailed investigations into the 
nonclassical properties of the microwave field in circuit 
QED experiments. 
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